Systematic study of Mn-doping trends in optical properties of (Ga,Mn)As 
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We report on a systematic study of optical properties of (Ga,Mn)As epilayers spanning the wide 
range of accessible substitutional Mnca dopings. The growth and post-growth annealing procedures 
were optimized for each nominal Mn doping in order to obtain films which are as close as possible 
to uniform uncompensated (Ga,Mn)As mixed crystals. We observe a broad maximum in the mid- 
infrared absorption spectra whose position exhibits a prevailing blue-shift for increasing Mn-doping. 
In the visible range, a peak in the magnetic circular dichroism blue shifts with increasing Mn- 
doping. These observed trends confirm that disorder-broadened valence band states provide a 
better one-particle representation for the electronic structure of high-doped (Ga,Mn) As with metallic 
conduction than an energy spectrum assuming the Fermi level pinned in a narrow impurity band. 

PACS numbers: 74.20.Mn, 74.25. Nf, 74.72.Bk, 74.76. Bz 



The discovery of ferromagnetism in (Ga,Mn)As above 
100 K [1 opened an attractive prospect for exploring the 
physics of magnetic phenomena in doped semiconduc- 
tors and for developing advanced concepts for spintron- 
ics. Assessment of a wide range of magnetic and trans- 
port properties of the material [2-4 showed that in ferro- 
magnetic (Ga,Mn)As with Mn dopings x > 1%, disorder- 
broadened and shifted host Bloch bands represent a use- 
ful one-particle basis for describing this mixed-crystal de- 
generate semiconductor. The common kinetic-exchange 
model implementation of this valence band theory and 
the more microscopic tight-binding Anderson model or 
ab-initio density functional theory can all be shown [5 to 
be mutually consistent on the level of atomic and orbital 
resolved band structure. The main utility of valence band 
theories have been in providing a qualitative and often 
semi-quantitative description of phenomena originating 
from the exchange split and spin-orbit coupled electronic 
structure and in assisting the development of prototype 
spintronic devices [4 . Other basic physical properties of 
(Ga,Mn)As, namely those reflecting the vicinity of the 
metal-insulator transition and localization and electron- 
electron interaction effects, remain to be fully understood 
and require to go beyond the commonly employed per- 
turbative or disorder averaged Bloch-band theories. 

In the insulator non-magnetic regime {x <C 1%), the 
system is readily described by localized Fermi level states 
residing inside a narrow impurity band separated from 
the valence band by an energy gap of magnitude close 
to the isolated MuQa impurity binding energy. Recently, 



a debate has been stirred by proposals, based in partic- 
ular on optical spectroscopy measurements [6 , that the 
narrow impurity band persists in high-doped (Ga,Mn)As 
with metallic conduction. Several phenomenological vari- 
ants of the impurity band model have been proposed 
for the high-doped regime [BHIO] which are mutually in- 
consistent from the perspective of the assumed atomic 
orbital nature of the impurity band states [5 . Fur- 
ther theoretical inconsistencies arise when recreating the 
phenomenological models microscopically with the con- 
straint of the experimentally determined moderate bind- 
ing energy of an isolated MuQa of ^a = 0-1 meV. A 
detached narrow impurity band does not persist in any 
of the variants of the model to dopings x > 1% corre- 
sponding to (Ga,Mn)As with metallic conduction [5]. 

The goal of this paper is to resolve the contro- 
versy between valence and impurity band approaches 
to (Ga,Mn)As from the experimental perspective by as- 
sessing doping trends in the optical spectra over a wide 
Mn concentration range in a consistently and control- 
lably prepared set of materials. We have optimized the 
growth and post-growth annealing procedures individu- 
ally for each nominal doping in order to minimize the 
density of compensating defects and other unintentional 
impurities and to achieve high uniformity of the epilay- 
ers in the growth and lateral directions, as detailed in 
the Supplementary material [TT. The trends observed 
in the unpolarized infrared and visible-range magneto- 
optical spectroscopies of this systematically prepared set 
of samples support the disordered valence band approach. 



They contradict previous conclusions which have been 
based on experiments in a more hmited set of (as-grown 
and annealed ) samples and which have been regarded as 
the key evidence for impurity band conduction in metal- 
lic (Ga,Mn)As [S [71 El EO] . 
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FIG. 1: (a) Infrared absorption of a series of optimized 
(Ga,Mn)As/GaAs epilayers with nominal Mn doping x — 
0.1 — 14% plotted from the measured optical transmissions 
of the samples (T) and of the reference bare GaAs substrate 
(To). Spectra of the 100 nm thick samples with x < 1% 
were divided by 5 for the consistency with those measured 
for the 20 nm thick higher doped samples, (b) Real part of 
the ac conductivity (lines) obtained from the measured com- 
plex conductivity in the terahertz range (points) and from 
fitting the complex conductivity in the infrared range to the 
measured transmissions, (c) Comparison of the infrared ab- 
sorption in as-grown and annealed 4.5% doped sample. Inset: 
Comparison to GaAs:C samples with carbon doping densi- 
ties 2 X 10^^ cm"^ and 2 x 10^° cm"^. (d) Height of the 
(Ga,Mn)As mid-infrared absorption peak as a function of Mn 
doping, (e) Position of the peak inferred from the trans- 
mission measurements and from the fitted ac conductivities. 
(f ) Zero frequency conductivities obtained from dc transport 
measurements and from extrapolated optical ac conductivi- 
ties measured in the terahertz range. 



Nominal dopings of our set of (Ga,Mn)As epilayers 
grown on GaAs substrates span a range from paramag- 
netic insulating materials with x < 1% to materials with 
X up to ^ 14%, corresponding to ^ 8% of uncompen- 
sated MuQa, and ferromagnetic transition temperatures 
reaching 190 K [11 . Samples with x < 1% have thickness 
of 100 nm. All epilayers with x > 1.5% are 20 nm thick. 



In these materials with large Mn doping, high quality 
epilayers are obtained only for thicknesses larger than 
~ 10 nm and lower than ~ 50 nm. All samples within the 
series have reproducible characteristics with the overall 
trend of increasing Curie temperature (in the ferromag- 
netic films), increasing hole concentration, and increasing 
magnetic moment density with increasing x. The sam- 
ples have a high degree of uniformity on a macroscopic 
scale as inferred from their sharp magnetic and transport 
singularities at the Curie point [11 . 

Samples used in transmission measurements have a 
polished back side of the wafers to minimize diffusive 
light scattering. The unpolarized transmission experi- 
ments ^t on the (Ga,Mn)As/GaAs samples and on the 
control bare GaAs substrate were performed at 300 K in 
the range 25-11000 cm~^ (3-1360 meV) using the Fourier 
transform infrared spectroscopy. The range between 120- 
600 cm~^ (15-75 meV) with strong phonon response in 
GaAs substrate is excluded from the data. Measurements 
of the complex conductivity in the low-frequency range 
8-80 cm-i (1-10 meV) of the (Ga,Mn)As/GaAs wafers 
(with measurement of the bare GaAs substrate as a refer- 
ence) were performed by means of terahertz time-domain 
transmission spectroscopy [11]. Magneto-optical experi- 
ments [11 at 15 K in the near infrared to visible range 
970-20160 cm~^ (1.2-2.5 eV) were performed primarily 
in the reflection geometry because of the small epilayer 
thickness of the optimized materials. Control magneto- 
optical measurements in transmission were done on a 
230 nm thick (Ga,Mn)As epilayer with the GaAs sub- 
strate removed after growth by wet etching. 

First we discuss the observed broad maximum in 
the unpolarized infrared absorption near 200 meV. In 
Fig. 1(a) we plot experimental data obtained directly 
from the measured infrared transmissivities. Real part 
of the ac conductivity curves shown in Fig. 1(b) rep- 
resent the best fit to the measured transmission in the 
THz and infrared ranges ^llj. The fit is anchored at 
low- frequencies by the directly measured THz conduc- 
tivity. (The scatter in the measured THz conductiv- 
ity reflects the precision of these measurements which 
is limited primarily by the quality of sample surfaces 
pT].) The position of the mid-infrared absorption peak 
in both representations of the measured data has a pre- 
vailing blue-shift tendency with increasing doping. This 
is reminiscent of the blue-shift of this spectral feature 
seen in our (and previously studied [T2]) control GaAsiC 
materials, shown in Fig. 1(c). We recall that for the non- 
magnetic hydrogenic acceptors it is established that the 
peak originates from transitions inside the semiconduc- 
tor valence band. Based on microscopic valence band 
theory calculations, these transitions have been also pre- 
dicted to yield the broad maximum in the mid-infrared 
spectra of (Ga,Mn)As which is rather insensitive to the 
exchange splitting of the valence band [I3l E]. Apart 
from the predicted prevailing blue shift, the peak po- 



sition in (Ga,Mn)As can have a non-monotonic depen- 
dence on doping as a consequence of momentum non- 
conserving transitions allowed by the strong disorder in 
(Ga,Mn)As [15]. 

A red shift of the measured mid-infrared peak reported 
in Ref. [7] in as-grown and annealed samples arranged 
by expected increasing hole concentrations has been pre- 
sented as the key evidence of the failure of the valence 
band theories. The data were phenomenologically inter- 
preted in terms of impurity band conduction persisting to 
the high-doped metallic (Ga,Mn)As. In the subsequent 
analysis it has been pointed out [15 , however, that the 
association of this peak to an impurity band is implausi- 
ble (i) because of the absence of the thermally activated 
dc-transport counterpart in the high-doped samples, (ii) 
because of the initial blue-shift of this mid-infrared fea- 
ture with respect to the impurity band transition peak 
in the very dilute insulating samples, and (iii) because 
of the appearance of the peak at frequencies above 2Ea 
which is the expected upper bound for impurity band 
transitions. Our data in Fig. 1 corroborate the conclu- 
sions of Ref. |T5j by demonstrating that the red-shift of 
the mid-infrared peak with increasing MuQa doping is 
not the general trend in (Ga,Mn)As materials prepared 
with the minimized number of compensating and other 
unintentional impurities. 

In Figs. 1(a) and (b) we have not included measure- 
ments in the as-grown epilayers because their characteris- 
tics are not fully reproducible, are more ambiguous, and 
the materials are less uniform [11 . Nevertheless, to make 
a connection to previous studies we have measured the 
mid-infrared peak in a nominally 4.5% Mn-doped sam- 
ple before and after annealing. As shown in Fig. 1(c), 
we observe a blue shift of the mid infrared peak after 
annealing, i.e., the red-shift is not observed in our ma- 
terials even if we use annealing to increase the effective 
doping. The same trend was confirmed in our as-grown 
and annealed samples with 12% nominal Mn doping. 

In Fig. l(d)-(f) we highlight the observed correlation 
between the peak position and its height, as well as the 
expected correlation from the valence band theories be- 
tween the amplitude of the peak and the dc conductivity. 
Note that in (Ga,Mn)As, the dc conductivity initially in- 
creases with doping because of the increasing number of 
itinerant holes provided by the Mn acceptors. At large 
dopings, the competition of increasing hole concentration 
and decreasing mobility yields a non-monotonic dc and 
ac conductivity trend (see also Ref. [11]). We also point 
out the close correspondence shown in Fig. 1(f) between 
the dc conductivities obtained from extrapolated THz 
data and from dc longitudinal and Hall transport mea- 
surements, which confirms the consistency of our optical 
data. 

Unlike the unpolarized optical spectroscopy, the 
magneto-optical effects are very sensitive to the mag- 
netic state of the system, their interpretation is more 
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FIG. 2: (a) Reflection MOD measurements in samples 
from the optimized series of 20 nm thick epilayers span- 
ning the whole studied range of Mn dopings of ferromagnetic 
(Ga,Mn)As. (b) Experimental position of the MCD peak 
(negative in the considered sign convention) as a function of 
nominal doping, (c) Comparison of the MCD measurements 
in the 12.5% doped as-grown and annealed epilayer. (d) The- 
oretical position of the MCD peak as a function of effective 
doping of uncompensated Mnoa impurities (see Ref. JJj); in- 
set shows theoretical MCD spectrum for 4% effective doping. 
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FIG. 3: Measured reflection MCD data in samples with dif- 
ferent (Ga,Mn)As epilayer thickness and nominal Mn doping 
1.5% (a), 12.5% (b). (c) Theoretical thickness dependence of 
the MCD spectra at 1% effective doping, (d) Transmission 
and reflection MCD measurements in a 7% nominally doped 
230 nm thick free standing (Ga,Mn)As epilayer. 



microscopically constrained, and they yield sharper spec- 
tral features. This is particularly valid for the infrared 
magneto-optical spectroscopy for which previous study 
[16] has shown that the disordered valence band the- 
ory with kinetic-exchange-split bands of ferromagnetic 



(Ga,Mn)As accounts semi-quant itatively for the overall 
characteristics of the measured data. We now extend the 
analysis to higher energies including transitions across 
the semiconductor band gap. 

The magnetic circular dichroism in the reflection geom- 
etry is given by MCD(deg) = 90/7r{R-^-R-)/{R-^^R-), 
where R'^^~^ is the reflected intensity of the light with 
the angular momentum +1 (-1). The main result of 
MCD measurements in our (Ga,Mn)As materials, shown 
in Figs. 2(a) and (b), is the observed large blue-shift of 
the MCD peak (the peak is negative in the above sign 
convention). Our data contradict the conclusion based 
on experiments in Ref. [9 which stated that the position 
of the MCD peak was independent of doping and that 
this was a signature of the pinning of the Fermi level in a 
rigid narrow impurity band [S] |10] . As in the case of the 
unpolarized infrared optical measurements, we have per- 
formed additional MCD experiments in pairs of as-grown 
and annealed (Ga,Mn)As materials. Upon annealing, we 
again observe a blue-shift of the MCD peak, as illustrated 
in Fig. 2(c). 

Theoretical modeling of the MCD spectra within the 
valence band theory, described in detail in Ref. fllj, is 
shown in Fig. 2(d). The magnitude of the peak and its 
absolute position depend on the detailed implementation 
of the model, e.g., on the way how disorder or band-gap 
renormalization effects are treated [11]. On a qualitative 
level and independent of these specific implementations 
of the model, the valence band theory with the antiferro- 
magnetic p — (i kinetic exchange [3 reproduces the overall 
experimental shape of the MCD spectrum, the sign of the 
MCD peak, and the blue shift of the peak with increasing 
doping. 

Because our experiments were performed in the reflec- 
tion geometry we have prepared several control samples 
to assess the role of multiple reflections. In Figs. 3(a),(b) 
we show a comparison of MCD spectra in (Ga,Mn) As epi- 
layers of different thicknesses controlled during growth 
or post-growth by etching [TT]. In the control thicker 
(Ga,Mn)As samples, the MCD peak is superimposed on 
oscillations caused by the multiple reflections. The os- 
cillations are completely suppressed in the 20 nm thick 
films and the reflection MCD peak position is therefore 
not affected by the film thickness in our series of opti- 
mized 20 nm thick (Ga,Mn)As epilayers. The transition 
from a thin film to a film where multiple reflections are 
important is captured qualitatively by the valence band 
theory calculations, as shown in Fig. 3(c). 

As an additional consistency check of the sign of our 
MCD signals with respect to previous MCD experiments 
in (Ga,Mn)As we have also measured the free-standing 
230 nm thick 7% Mn-doped (Ga,Mn)As epilayer in both 
reflection and transmission MCD geometry. In the latter 
case the dichroism is given by MCD(deg) = 90 /tt (T+ — 
T-)/{T^^T-), where T+(-) is the transmitted intensity 
of the light with the angular momentum -fl (-1). Our 



transmission MCD data shown in Fig. 3(d) are very sim- 
ilar to measurements in samples with comparable doping 
and film thickness reported in Ref. pT] , 

To conclude, we have not observed the signatures in 
the optical spectra which in previous studies have been 
interpreted as general trends in metallic (Ga,Mn)As and 
as an evidence of Fermi level states residing in a nar- 
row impurity band detached from the valence band. Our 
unpolarized infrared and visible-range magneto-optical 
measurements performed on an extensive series of sys- 
tematically prepared (Ga,Mn)As materials support the 
valence band theory representation of the effective one- 
particle spectrum of high-doped (Ga,Mn)As with metal- 
lic conduction. 
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SUPPLEMENTARY MATERIAL 



INTRODUCTION 

The supplementary material describes a series of (Ga,Mn)As epilayers intended to provide a basis for systematic 
experimental studies of this ferromagnetic semiconductor. For each nominal Mn doping x, the growth and post-growth 
annealing conditions were separately optimized in order to achieve the highest Curie temperature Tc attainable at 
the particular x. The highest Tc criterion was found to lead simultaneously to layers with maximized uniformity and 
minimized compensation by unintentional impurities and defects. 

The nominal doping x spans a wide range from non- ferromagnetic insulating materials with x < 1% to the highest Tc 
materials with x ~ 13%. Within this nominal Mn doping range and for film thicknesses smaller than approximately 
50 nm and larger than 10 nm in case of higher doped samples with x > 1.5, we obtained a series of (Ga,Mn)As 
materials with reproducible characteristics. They show an overall trend of increasing saturation moment (ascribed to 
substitutional MuQa) with increasing x , increasing T^ and increasing hole density. The materials have no apparent 
charge or moment compensation of the substitutional MuQa impurities and have a large degree of uniformity reflected 
by sharp magnetic and transport singularities at Tc. 

Because experimentally it is difficult to determine on a precise quantitative level all relevant parameters of a given 
material and also because accurate quantitative theoretical modeling of even an idealized (Ga,Mn)As is not feasible 
we are convinced that a meaningful examination of intrinsic properties of (Ga,Mn)As ferromagnetic semiconductor 
materials should be performed on the level of qualitative or at most semi-quantitative trends observed over a wide 
doping range in a consistent set of materials. Our series of epilayers represents an attempt to provide such an 
experimental basis. The optimization of the materials in the series, which is performed individually for each nominal 
doping, minimizes the uncertainties in the experimental sample parameters and produces high quality epilayers which 
are as close as possible to to uniform uncompensated (Ga,Mn)As mixed crystals. We prefer this method for studying 
the doping trends over, e.g., comparisons of as-grown and annealed samples because the as-grown state of a sample 
is not well defined; the characteristics of as-grown samples are not well reproducible, the materials are less uniform 
and can have a large number of unintentional interstitial Mn impurities which can compensate both the local moment 
and holes produced by MuQa and whose density is difficult to measure accurately. Similarly, the approach based on 
individual optimization at each nominal doping provides a qualitatively better control over the whole set of sample 
parameters and therefore a better approximation to idealized (Ga,Mn)As materials than the combinatorial method 
which on one physical wafer produces a wide range of samples due to the stoichiometry gradients across the wafer 
imposed during a single growth experiment. For further discussion highlighting the importance of the individual 
optimization of material synthesis at each nominal doping see also Ref. pj . 

GROWTH 

The growths were done in a Veeco Gen II molecular-beam epitaxy (MBE) system. We used 2-inch semi-insulating 
0.5 mm thick GaAs substrates mounted on a molybdenum In- free sample holder by mechanically fixing the sample 
at the edges. Each growth started by a standard high-temperature grown GaAs buffer of about 200 nm thickness. 
After the growth interruption the substrate was cooled down for the low-temperature growth of (Ga,Mn)As; during 
the cooling period the As source was closed when the substrate temperature dropped below 400° C. 

The (Ga,Mn)As layers were grown at the growth rate of 0.2 monolayers/second, corresponding to temperatures of 
940° and 920° C at the tip and the base zone, respectively, of the dual filament Ga source. The Mn flux, and thence 
X, was determined by measuring the ratio of the beam equivalent pressures (BEP) of Mn and Ga sources before each 
growth. The Mn content was cross-checked on several samples by secondary ion mass spectroscopy (SIMS) and by 



comparing the growth rates of GaAs and (Ga,Mn)As; the results of the three methods agreed within Ax ~ 0.015. The 
BEP-based method yields the best reproducibility and guarantees the correct ordering by increasing x even though 
the relative uncertainty in x may be large especially at low dopings. 

The Arsenic was supplied by a valved cracker cell in the As4 regime (cracking zone temperature 600° C). The As- 
to-(Ga+Mn) flux ratio was calibrated using the As-controlled RHEED intensity oscillations. For this purpose the 
substrate temperature was first set to 500° C. Next, the As source was closed and the Ga shutter was simultaneously 
opened for 30 - 40 s allowing accumulation of several Ga monolayers on the substrate surface. After closing the Ga 
shutter and re-opening the As source several growth oscillations were observed in the RHEED image with period 
corresponding to the As flux. We have checked that this technique yields essentially the same calibration as the 
As-rich to Ga-rich surface reconstruction transition method used in Ref. 1 . 




260°C: polycryst. 240°C: 3D 220°C: 2D 

FIG. 4: RHEED images taken during growths of (Ga,Mn)As materials with nominal 7% Mn doping, stoichiometric flux ratio 
and different substrate temperatures Tg. The 2D growth mode is achieved only at sufficiently low Tg. 

An extremely critical growth parameter of (Ga,Mn)As is the substrate temperature. For its measurement we used 
the GaAs band edge spectrometer (kSA BandiT) mounted on the central pyrometer port, normal to the substrate. 
In this arrangement, the thermal radiation of the heater has sufficient intensity near the band gap wavelength to 
serve as a radiation source for the transmission measurement even at temperatures below 200° C, provided that no 
diffuser plate is inserted between the substrate and the heater. In Fig. [4] we show examples of RHEED patterns for 
the nominally x = 7% doped (Ga,Mn)As material grown at stoichiometric 1:1 ratio of As:(Ga+Mn) and for different 
substrate temperatures Tg. For this nominal doping and flux ratio, the transition between 3D and 2D growth modes 
occurs for Tg between 240 and 220°C. In Fig. [s] we plot the 3D/2D growth boundary as a function of nominal Mn 
doping, reproducing essentially the previously reported observation [2] that the maximum allowed Tg for the 2D 
growth rapidly decreases with increasing doping. It can also be seen in Fig. p] that the 2D/3D transition temperature 
depends on the As:(Ga+Mn) ratio. The excess As flux compensates for the higher As re-evaporation rate and allows 
thus to increase the growth temperature without initiating the 3D growth mode. On the other hand, this gain in the 
maximum growth temperature, positive from the perspective of the epitaxial growth, does not fully compensate the 
increase in the density of AsQa antisites, and the highest quality (Ga,Mn)As layers are still obtained when growing 
(i) in the 2D growth mode, (ii) as close as possible to the 2D/3D transition line, (iii) with the stoichiometric ratio 
As:(Ga+Mn)^l:l. 

FILM THICKNESS AND POST-GROWTH ANNEALING 

Film thickness and post-growth annealing conditions are other important factors determining the quality of the 
resulting (Ga,Mn) As materials. In Fig. |6] we show the dependence of the Curie temperature Tc on the annealing 
time for three different annealing temperatures for our record Tc = 188 K sample with nominal 13% Mn doping and 
film thickness 15 nm. These curves illustrate the common trend in annealing (at temperatures close to the growth 
temperature) suggesting two competing mechanisms. One yielding the increase of Tc and ascribed in a number of 
reports to the removal of charge and moment compensating interstitial Mn impurities (see also our detailed annealing 



As:(Ga+Mn)=3:1 

As:(Ga+Mn)=2:1 

As:(Ga+Mn)=1:1 




X (%) 

FIG. 5: Boundary in substrate growth temperature separating the 3D and 2D growth modes as a function of the nominal doping. 
Bottom (black), middle (blue) and upper (red) curves correspond to As:(Ga+Mn) ratio of 1:1, 2:1, and 3:1, respectively. Empty 
circles denote samples which retained the 2D growth mode at thickness 50 nm and more, full diamonds denote samples which 
became 3D at thicknesses between 20 and 50 nm. 
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FIG. 6: (a) Curie temperature as a function of annealing time for three different annealing temperatures for a 15 nm thick 
(Ga,Mn)As epilayer with nominal 13% Mn doping; the record Tc — 188 K was achieved using annealing in air at 160° C. (b) 
Conductivity (upper curve at 4.2 K and bottom curve at room temperature) monitored simultaneously during annealing at 
160° C of the same epilayer. 



studies in Ref. [3]). In highly doped samples [x 



> 



a second mechanism is observed which starts to slowly 



reduce Tc after sufficiently long annealing times and at higher annealing temperatures. The origin of this detrimental 
mechanism is not established but can be ascribed to a partial removal of substitutional Mn or to a formation of 
some stable complexes (e.g. 2MnGa+Mnint clusters) in the film. The lower annealing temperature the longer time 
is required for observing the downturn of Tc but also longer time is required to reach the maximum Tc. Because of 
the two competing mechanisms, the absolutely highest Curie temperature for the given nominal doping is achieved 
at some intermediate annealing temperature and time, as illustrated in Fig. |6Fa). 

We also note that while the increase of Tc with annealing is closely correlated with the increase of the conductivity 
of the film, which is often used as a practical tool for monitoring the annealing progress, the maximum Tc is typically 
reached at a later annealing time than the maximum of the conductivity, as illustrated in Fig. [6Fb). This effect is 
stronger in very thin samples and can be ascribed to layer thinning due to oxidation during the annealing. The 
difference between the peak conductivity and the conductivity at maximum Tc remain, however, small at practical 
annealing times. 

For a given nominal doping, the highest attainable Tc is reached only in thin films, typically thinner than ~50 nm. 
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FIG. 7: (a) Annealing study of a 100 nm grown 13% doped epilayer directly annealed (black curve), annealed after thinning 
(red curve) and partially annealed, thinned and re-annealed (blue curve), (b) Resistivities p{T) and the dp/dT singularities of 
50 and 20 nm films with the same nominal Mn doping of 7% showing good quality of optimally grown and annealed 50 nm 
films but not reaching fully the quality of optimized thinner films. 
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FIG. 8: Upper panel shows temperature dependent magnetization M(T) and dp/dT illustrating high magnetic quality of the 
optimized films; resistivity p(T) is shown in the lower panel. Singularities in M{T) and dp/dT can be used to infer the Critical 
temperature while the broad maximum in p{T) is above Tc. This is a common result illustrated in this figure on an optimized 
11% nominally doped (Ga,Mn)As film. 



(Note that in ultra-thin films with thicknesses below ~10 nm, which start to show a clear enhancement of resistivity, 
Tc also does not reach the highest values obtained at the given nominal doping in the thicker films.) In Fig. ITFa) we 
illustrate the importance of the film thickness for obtaining the high quality (Ga,Mn)As materials. A 100 nm thick 
film is grown with nominal 13% doping and, unlike the thin record Tc film discussed above, here the maximum Tc 
achieved by annealing is only about 140 K. However, even after passing this maximum the Tc can still be increased 
when the annealed 100 nm thick film is thinned down (to, e.g., 20 nm) by wet etching and annealed further, as also 
shown in Fig. JTFa). The highest Tc is achieved when the film is thinned down before any annealing. 

Our systematic studies of (Ga,Mn)As materials therefore focus on films with thickness below 50 nm and the most 
extensive series of carefully optimized materials was prepared using 20 nm thick epilayers. To illustrate that the 
50 nm films can have good quality but do not fully reach the Tc^ low resistivity, and sample uniformity after optimal 



annealing as the 20 nm films with the same nominal doping we show in Fig. iT^b) a comparison of resistivities, p(T), 
and the dp/dT singularities at the critical point for the 50 and 20 nm thick films (with 7% nominal doping). The 
resistivity is smaller and the dp/dT singularity is sharper in the 20 nm film confirming a higher degree of uniformity. 
The position of the dp/dT singularity, corresponding to Tc (see Fig. [sland Ref. [4 for details), is shifted to slightly 
higher temperature in the thinner film. 

SERIES OF (Ga, Mn) As MATERIALS PREPARED UNDER OPTIMIZED GROWTH AND ANNEALING 

CONDITIONS 
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FIG. 9: Measurements of the insulator to metal transition induced by increasing Mn doping in (Ga,Mn)As. (a) Temperature 
dependent conductivities in the insulating low doped and metallic high doped samples, (b) Red curve: conductivity at 4.2 K 
as a function of doping showing abrupt increase of the conductivity at x ~ 1.5% and only a weak dependence (on this large 
scale) of conductivity on doping beyond this value, (b) Black curve: dependence of Curie temperature on doping. 

After finding the optimal growth and post-growth conditions for each individual nominal doping we obtained a 
series of samples whose characteristics are summarized in Tab. [Tj The samples can be divided into several groups: at 
nominal dopings below ~ 0.1% the (Ga,Mn)As materials are paramagnetic, strongly insulating, showing signatures 
of the activated transport corresponding to valence band - impurity band transitions at intermediate temperatures, 
and valence band - conduction band transitions at high temperatures (see upper panel in Fig. [9]). For higher nominal 
dopings, 0.5 ^ x ^ 1.5%, no clear signatures of activation from the valence band to the impurity band are seen 
in the dc transport, suggesting that the bands start to overlap and mix, yet the materials remain insulating. At 
X ^ 1.5%, low-temperature conductivity of the film increases abruptly by several orders of magnitude (see bottom 
panel of Fig. [9|, and the system turns into a degenerate semiconductor. The onset of ferromagnetism occurs already 
on the insulating side of the transition at x ~ 1% and the Curie temperature then steadily increases with increasing 
doping up to nominal doping of ~ 13%. Beyond these dopings we did not succeed in introducing more substitutional 
MuGa moments into the system. 



In Fig. 10 we show the dp/dT measurements throughout the series of optimized samples spanning the nominal 
doping range from ^1.5 to 13% which illustrate the high quality of all the epilayers within the series. 

Fig. [TT] shows results of the Hall effect measurements at 4.2 K. For this purpose the samples were lithographically 
patterned into Hall-bars of 60/im width. It can be seen in Fig. 11 that the Hall signal is affected by longitudinal 



magnetoresistance of the samples, especially at low dopings. Therefore, we extracted p from high field data and by 
fitting the measured transversal resistance Rxy by 



Rxy = B/{epd) + kiRj^a, + k2Rl 
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FIG. 10: dp/dT measurements illustrating high magnetic quality of the optimized films throughout the whole series spanning 
the nominal doping range from ^1.5 to 13% and Tc's from ^25 to 186 K. In the inset detail of the dp/dT cusp with maximum 
normalized to unity is shown. 
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FIG. 11: Transport measurements of the series of optimized materials, (a) Longitudinal resistance of selected samples (nor- 
malized to Rxx at i? = 0). (b) Hall resistance; note the strong anomalous Hall-effect contribution. 



where d is the sample thickness and ki and /c2 are fitting constants reflecting the anomalous Hall effect and possible 
imperfections in the geometry of the Hall bars. 

We also emphasize that, apart from the common experimental scatter and from the corrections due to the non- 
zero magnetoresistance and due to the anomalous Hall effect, the carrier density can in principle be inferred only 
approximately from the slope of the Hall curve in a multi-band, spin-orbit coupled exchange-split system such as the 
(Ga,Mn)As. The error bar due to the multi-band nature is estimated to ~ 20% [5 . Due to these uncertainties we can 
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FIG. 12: (a) Hole density at 4.2 K extracted from the high-field slope of the Hall resistance, (b) Hole mobility at 4.2 K. Inset: 
Hole mobility vs. hole density. 




FIG. 13: Magnetization measurement of the series of optimized materials. Saturated moment is plotted by green squares and 
refers to the left y-axis. Right y-axis translates the saturated moment into the density Nmu of uncompensated substitutional 
Mn atoms. Black circles denote the hole density p inferred from the Hall effect. Dashed line indicates the case when the 
effective uncompensated doping equals to the nominal doping x, i.e. when Nmu — p — 4x/a|^. 



only make semi- quantitative conclusions based on the measured Hall effect hole densities. Among those conclusions 
is the confirmed overall trend of increasing p (and decreasing mobility fi) with increasing doping in the optimized 



materials, as shown in Fig. 12 



Magnetization measurements are summarized in Fig. [13] The figure confirms the increasing trend in the saturation 
moment (corresponding to the density of uncompensated MuQa moments) with increasing nominal doping up to 
X ^ 13. Assuming 4.5 fiB per Mn atom [6 we can estimate the density Nmu of uncompesated MuQa from the 
magnetization data, as shown by the right y-axis in Fig. 13 An important conclusion can be drawn when comparing 



this estimate with the hole density p previously estimated from the Hall resistance. Since there is no apparent deficit of 
p compared to Nmui ^ind since the interstitial Mn impurity compensates one local moment but two holes we conclude 
that interstitial Mn is completely (within the experimental scatter) removed in our optimally annealed epilayers. The 
series of optimized materials is essentially free of charge and moment compensation and the density of both holes and 
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local moments is given by the density of Mn in substitutional sites. As can be inferred from Fig. [13J the maximum 
effective doping of substitutional Mn atoms reaches ^ 8% at nominal doping x ^ 13%. 
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1.03 
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20 
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- 
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F056 
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20 
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- 
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TABLE I: Table summarizing the basic characteristics of samples from the series of optimized materials. 



OPTICAL EXPERIMENTS 

The infrared transmission spectra were recorded at 300 K using beam condenser kit on a Nicolet 6700 FTIR 
spectrometer with 4 cm~^ resolution and Happ-Genzel apodization. The spectra in the MIR region (350-4000 cm~^) 
were recorded using KBr beamsplitter and DTGS detector; the spectra in the NIR region (4000-11000 cm~^) were 
recorded using CaF2 beamsplitter and PbSe detector. The spectra in the FIR region (25-700 cm~^) were recorded on 
an evacuated Bruker IFS 113v spectrometer with 2 cm~^ resolution and Blackman-Harris 3- Term apodization using 
6 /im thick mylar beamsplitter and a DTGS/PE detector. 

For the THz time-domain experiments we used the experimental setup described in Ref. [7 . The linearly polarized 
broadband THz pulses formed a ~ 3 mm spot on the sample and probed the in-plane complex conductivity a of 
(Ga,Mn)As films. Each experiment consisted of two consecutive measurements: one of a time-domain signal Es{t) 
transmitted through the (Ga,Mn)As/GaAs wafer and another one of a reference signal Er{t) transmitted through a 
bare GaAs substrate. The Fourier transforms of the signals allow us to calculate the complex transmittance spectrum 
of the film: t{uj) = Es{uj)/Er{oj). The spectra of complex refractive index N (and of the conductivity a = —iuje^N'^) 
are then retrieved following the procedure described in Ref. [7 . In our samples with the substrate polished from the 
back side, the wafer thickness precision is 1-2 /im. This leads to the experimental error of the order of 30-60 r^~^cm~^ 
in the magnitude of the measured Re cr in the THz range which is plotted as dots in Fig. 1(b) in the main text. 

As emphasized in Ref. [8 , no general practical sum rule exists for the low- frequency part of the ac conductivity of 
(Ga,Mn)As films. Consequently, we did not attempt to interpret the shape of the observed low-frequency conductivity. 
We used a non-Drude (or generalized Drude) formula [9 which allowed us to obtain a very good fit of the data 
simultaneously in the THz and far-infrared frequency ranges. The conductivity peak in the mid-infrared region was 
fitted by adding a sum of three damped harmonic oscillators (in order to take into account its asymmetry) . The entire 
fitting formula then reads. 



(j{uj) 



O'dc 



[1 - {iujrY 



-IUJCq 



V^ fi 



zcjeoCo 



(2) 
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We obtained excellent simultaneous fits of the THz, far- and mid-infrared transmittance data by using the above 
formula (with ol ^ \ and /3 ^ 1) and the transfer matrix formalism ^10^ to calculate the transmission through 
(Ga,Mn)As/GaAs system. We emphasize that the uncertainty in the resulting mid-infrared peak position of Re a is 
much smaller than the shift of the peak observed across the studied Mn doping range. 

The spectra of MCD were studied by the experimental procedure described in [11 . In particular, the photoelastic 
modulator (PEM) was used to alter the polarization state of the light. The sign of MCD was determined from an 
independent experiment, where the helicity of the light was controlled by a quarter wave plate. All MCD experiments 
were performed at temperature of 15 K. The angle of incidence of the light beam was smaller than 1 deg (i.e., very close 
to the normal incidence). An external magnetic field of 540 mT was applied perpendicular to the sample plane during 
the experiment. The magnetic field direction-independent effects were removed by taking the difference between the 
data measured with +540 mT and -540 mT. The data were also corrected for the MCD signal from the cryostat 
windows. 

PARABOLIC MODEL OF OPTICAL DIELECTRIC TENSOR IN METALLIC GaMnAs 

The interband conductivity of a system can be calculated through the Kubo formula : 

where p^^ = l\a'k\pa\h'k) and £^5^ = ^hk — ^ak- Its relation to the dielectric tensor is 

UJ 

which obey the Kramers-Kronig relations 

im[..,] = -'-pr 5£Md.' = -^-P r ^-^^d.' (6) 

IT J_^ Uj' -UJ TT Jo ^' - ^^ 

Ignoring the dc-contributions to the dielectric constant, the dielectric tensor from inter-band contributions can be 
written as 

inter/ X , r ^TTe^fl^ ^ (/a,/c " fh.k) PgbPba (^^ 

This expressions amounts to saying that one requires e^^y^^{uj = 0) = and that ^^^^^{uj = 0) = 
^^^2y Xlk ai-h J^^^^ahPha Yii calculating the dielectric tensor it is simpler to calculate first the contributions 

containing delta functions and then obtain the other contributions through the KK relations. More specifically 

Im[4r(^)] = ^5^ E ^^''^''J'^'^ MPl.PUK^-Eba) (8) 



\c,ay^b 



and 



Re[4r(-)] = -^5^ E ^\2^'^'^ MP>L]^{n.^-Eba) (9) 

Kinetic- exchange parabolic band model The model that we consider is one of parabolic bands with effective masses 
'^ci'^hhi'^ihi'^so- The conduction bands are unoccupied and the valence bands are partially occupied with a given 



Fermi energy Ep measured from the top of the valence band without exchange splitting as shown in Fig. 14 Here the 



label c G 6 and hh^ Ih^ so G a. The exchange splitting (see Fig. 14) is given by A^^ = ^sd^b + ^pd^a- Here /S^sd is due 
to the ferromagnetic direct-exchange coupling between the conduction band 5-orbitals and the Mn (i-orbitals from 
the same Ga-sublattice of the host zinc-blende crystal; Ap^^ is due to the antiferromagnetic kinetic-exchange coupling 
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FIG. 14: Sketch of parabolic model (only two valence bands shown). 



which originates from the hybridization between the p-orbitals from the As sublattice and the Mn d-orbitals. Agd 
and Apd are defined as positive numbers. 

In what follows we consider pairs of valence bands and conduction bands (i.e. the transitions of the ±3/2/i/i to the 
±l/2c, etc.) with the labeling: 
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AsdSb + Eg 



Eba — Ea — Eb — — h h Ah 



2mr 



2m„ 



2iiab 



(10) 

(11) 
(12) 



where fiah = '^a'^cK'^a + ^c)- We will also use the usual notation for the dipole matrix elements Ep = 2moP^/fi^ 
andP = -(zn/mo)(5|p^|X). 

Then for the Imaginary component of the diagonal dielectric tensor the contribution from the transition from a to 
b ih = 0) is: 



lm[e,, (a;)J - ^^^^ r,,, y^^ 



3/2 
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2fiab 
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(tUOy 2Hab 



(13) 
(14) 



where we have used Ry = h^ /2mQa^ = e^/2ao and e^ /Aha^ = e^/87rfiao = 1.83 x 10^1^ "'^cm ^ . 
The real part of Cxx can then be found from the KK relations: 
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(16) 
(17) 

(18) 



Effects of disorder broadening. Ferromagnetic (Ga,Mn)As is a very highly doped and therefore strongly disordered 
semiconductor material. This affects the magneto-optical effects by introducing a finite life time broadening to 
the quasiparticles and by allowing momentum non-conserving transitions. The former can be easily approximated 
by introducing an imaginary part to the energy denominators, i.e., iS -^ iV (we used 100 meV). On the other 
hand, within our approximation, the effects of non-conserving momentum transitions translate in a lower bound of 
the transitions which we approximate by increasing the conduction band effective mass, rric ~ 0.7mo. Also, the 
band gap renormalization from the electron-electron interactions, considering the larger effective mass, is given by 
54 X (p[10^^cm~^])^/^ meV [12]. This leads to the above expressions being replaced by 
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Im[6Lr(^)]dis = -If^ d(W)Re[6Lr (/^c^O]^^^^" {ful^-f^y'l r^ ^^^) 
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